Quantum switches and quantum memories for matter-wave lattice solitons 
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We study the possibility of implementing a quantum switch and a quantum memory for matter 
wave lattice solitons by making them interact with "effective" potentials (barrier/well) corresponding 
to defects of the optical lattice. In the case of interaction with an "effective" potential barrier, 
the bright lattice soliton experiences an abrupt transition from complete transmission to complete 
reflection (quantum switch) for a critical height of the barrier. The trapping of the soliton in an 
"effective" potential well and its release on demand, without loses, shows the feasibility of using 
the system as a quantum memory. The inclusion of defects as a way of controlling the interactions 
between two solitons is also reported. 



I. INTRODUCTION 

Bose Einstein condensates (BEC) in optical lattices 
have attracted during the last years a lot of attention 
both in the mean field regime [l[ as well as in the strongly 
correlated regime Q. One of the main reasons for this 
activity burst is the high level of control achieved in the 
experiments of ultracold gases in optical lattices which 
permits to explore a broad range of fundamental phe- 
nomena. 

In the mean field regime, a huge interest has been de- 
voted to nonlinear dynamics of matter-waves in periodic 
media and specifically to matter-wave solitons. Matter- 
wave solitons are self stabilized coherent atomic struc- 
tures that appear in nonlinear systems due to the bal- 
ance between the nonlincarity and the dispersive effects. 
The nature of the solitons supported by Bose Einstein 
condensates (BEC) is determined by the character of 
the interactions: attractive (repulsive) nonlinearity sup- 
ports bright Q (dark Q) solitons. In the presence of 
an optical lattice, this scenario changes completely due 
to the appearance of a band structure in the spectrum 
and the possibility of having either bright or dark lattice 
solitons with either repulsive or attractive interactions 
arises. Very recently, the first experimental demonstra- 
tion of bright lattice solitons in repulsive condensates was 
reported @. 

Since the first proposals of BEC lattice solitons Q, 
there has been an explosion of contributions regarding 
generation, mobility and interactions of this novel type 
of matter-wave solitons both in one-dimensional systems 
0, H, [§] and in higher dimensions [l(|. The interest is 
mainly centered in bright matter-wave lattice solitons 
due to their potential applications in energy and informa- 
tion transport in nonlinear systems. The fact that mat- 
ter wave solitons are massive permits to generate them 
at rest and to move them after an appropriate transfer 
of momentum. Proposals for controlling the dynamics of 
bright gap solitons are mainly devoted to the manipu- 
lation of the optical lattice [ll| and to the modification 



of the nonlinearity [12|. Nevertheless, a complete control 
on the dynamics of bright matter-wave gap solitons also 
requires a profound knowledge of their interactions with 
defects. 



The interaction of solitons with local inhomogeneities 
is a subject that appears in the literature in different con- 
texts and has been studied in the framework of different 
nonlinear equations (see for instance [13|). In particu- 
lar, the nonlinear Schrodinger equation with point-like 
defects either in the continuum regime [l4| or in discrete 
systems [l5| has deserved special attention. Extended de- 
fects have also been addressed in this framework [lj| . In 
nonlinear optics, the coupled mode equations have been 
used to study collisions of moving Bragg solitons with 
finite size 

and point like defects Very 
recently, interactions with defects in the context of con- 
tinuous matter-wave solitons have also been addressed 

MM. 



In this paper we focus on bright lattice matter-wave 
solitons and propose different possibilities of control of 
their dynamics by making them interact with defects of 
arbitrary amplitude and width. Specifically we will show 
how to change the direction of movement (a complete 
bounce back) of the soliton and how it can be stored and 
retrieve on demand. In Sec. [Til the physical system con- 
sidered and the model used is introduced. In Section IIIII 
we present the results concerning the interaction of soli- 
tons with an "effective" potential barrier, where the pos- 
sibility of implementing a quantum switch arises. Next, 
in Sec. IIVI the results regarding the interaction with an 
"effective" potential well will be shown and the use of 
the system as a quantum memory will be discussed. The 
possibility of controlling the interactions between two lat- 
tice solitons by placing a defect at the interaction point 
is discussed in Sec. IVl We will conclude in Sec. I VII 
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II. PHYSICAL SYSTEM 



We consider a zero temperature 87 Rb condensate con- 
fined in a one-dimensional geometry and in the presence 
of an optical lattice. The description of the system is 
performed within the one dimensional Gross Pitaevskii 
equation (GPE): 

tftfJM = (_ A + V (x) + gty (x, t) | 2 ) V t) , 

(1) 

where g = 2ha s u) t , being a s the s-wave scattering length 
and tot the radial angular trapping frequency, is the av- 
eraged one-dimensional coupling constant. The external 
potential is given by: 



V(x) 



2 x 2 + y sin 2 (^), 



(2) 



which describes both the axial trapping potential, with 
angular frequency lu x , and the optical lattice, with spatial 
period d = A/2sin(0/2), being A the wavelength of the 
lasers forming the optical lattice and 8 the angle between 
them. The depth of the optical lattice, Vq, is given in 
units of the lattice recoil energy E r = h 2 k 2 /2m where 
k = ir/d is the lattice recoil momentum. 

The generation of the bright lattice soliton is per- 
formed as it is reported in [8j using parameters close 
to the experimental realizations @. The procedure is 
briefly summarized in what follows. The starting point 
is a 87 Rb condensate (a s — 5.8nm, m = 1.45 x 10~ 25 Kg) 
of N — 500 atoms in the presence of a magnetic trap with 
frequencies u>t — 715 x 2tt Hz and uj x = 14 x 2ir Hz and an 
optical lattice, with potential depth Vq = lE r and period 
d = 397. 5nm. The axial magnetic trap is suddenly turned 
off and the appropriate phase imprinting, corresponding 
to phase jumps of tt in adjacent sites, is performed @. 
After the phase imprinting, the system evolves towards a 
negative mass, self-maintained staggered soliton at rest 
centered at x — 0, which contains approximately 35% of 
the initial atoms (N = 187) and extends circa 11 sites. 
The exceeding atoms are lost by radiation. 

The total energy of the generated bright lattice soliton 
can be calculated using the energy functional of the GPE 
([T|) that contains the total kinetic (-E/T), interaction (Ei) 
and potential (E p ) energies: 




1 + e- fe ' 2 " 2 cos(2fc'a;o) , /2 
1- k 



2 V 



3 + e ~ 2fe ' 2,)2 cos(4fc'a;o) + 4e~^ cos(2k'x ) 



1 + e- k ' 2 ^ cos{2k'x ) - e~ k2r i 2 cos(2kx ) 



-e k - r ' 2 cos(2/c_x ) - ie fc +'' 2 cos(2k + x ) 



(4) 



where B = |A| 2 0Fr//4, k' = 2tt/X' and fc ± = k' ± k. 
Fixing N — 187, S — 0.07/xm and xq = we obtain 
a minimum of Eq. ([¥]) corresponding to 1.31-E r . Ex- 
act numerical integration of Eq. ^ gives a total en- 
ergy of 1.35E r . This shows the high level of accuracy 
that the used variational method provides. Also, an 
strong agreement is found when we evaluate each term 
of Eq.© by using the variational method [§| (exact in- 
tegration): El = 0.85E r (0.85E r ), E, = 0.12^(0.13^.) 
and E p = 0.34^(0.37^). 

By calculating the linear band spectrum of the system 
we obtain an energy at the band edge of 1.25i? r . This 
value is in good agreement with the total energy obtained 
previously (Eq.Q) without the non linear term. The 
linear band spectrum predicts an effective mass at the 
edge of the first Brillouin zone corresponding to m e / / = 
—0.15m. 

Once the lattice soliton is created, it is set into mo- 
tion by applying an instantaneous transfer of momen- 
tum at t = 0. It has to be large enough to overcome 
the Peierls-Nabarro (PN) barrier @, H2] but sufficiently 
small to assure that the soliton remains in the region of 
negative effective mass, i.e., 0.009/c/i < p < 0.2kh. The 
soliton starts to move opposite to the direction of the 
kick, manifesting thus its negative effective mass. The 
concept of the effective mass is used through the paper 
to give an intuitive explanation of the observed dynam- 
ics. Nevertheless, all the results presented in what follows 
have been obtained, without any approximation, by di- 
rect integration of Eq. |1| . 

At a certain distance x rn of the initial position of the 
soliton (x = 0), the lattice potential V(x) is modified in 
the following way: 



E 



fi^)i 4 



V(x)M x)\'< 



dx. 



As clearly observed in the numerical solutions of the 
GPE, the density profile of the bright lattice soliton 
inside each well is shifted with respect to the mini- 
mum of the optical lattice. This shift has to be taken 
into account to properly calculate the energy. Assum- 
ing a constant shift 5, the Ansatz ip(x) — Aexp(— (x — 
x ) 2 /2r] 2 ) cos(27rx/(A')), with A' = 2d + 5, A the ampli- 
tude and r\ the width of the Gaussian envelope [9j , leads 
to the following energy functional: 



f Vq sin 2 (2f ) + V m (l - {X 2 X J )2 ) if x m -l/2<x<x m + 1/2; 
1 Vbsin 2 (^) otherwise. 

where V m can be either positive or negative and a = 6d 
for all the cases. For V m < 0, the local decrease of the 
lattice potential corresponds to an "effective" potential 
barrier for the soliton due to its negative effective mass 
while if V m > 0, i.e. a local increase of the potential 
acts as an "effective" well for the soliton. For the case 
of an effective barrier, x m is fixed to match exactly with 
a minimum of the optical lattice while in the effective 



3 



well case, x m corresponds to a maximum of the optical 
lattice. We have checked that the results reported in the 
following do not depend strongly on the specific shape 
of the potential by reproducing them with Gaussian and 
square potentials for the defect. 

To analyse the interaction of the bright lattice soliton 
with the defect, it is crucial to know the total energy 
of the soliton while it moves. The applied variational 
Ansatz with the shift in the periodicity is meaningful only 
in the static case and cannot be used to study soliton dy- 
namics, since the soliton is always chirped with respect 
its center Q. Therefore, to study dynamical behaviour 
one has to rely on numerical simulations. We have nu- 
merically calculated the contributions to the total energy 
of the soliton as a function of time. Immediately after 
the kick, the soliton expels atoms and its energy abruptly 
decreases becoming much smaller than the energy that it 
would need to remain at rest at the edge of the first Bril- 
louin zone. In the framework of the linear band theory 
this would correspond to displace the particle from the 
edge of the band of the first Brillouin zone by changing 
its quasimomentum. To illustrate the dynamics of the 
system, we consider the case in which we give a kick of 
p = O.lkh to the generated soliton at rest. At t = 0, just 
after the kick, the total energy of the soliton is its energy 
at rest plus the contribution of the transfer of momen- 
tum, i.e. E = \:ibE r + (0.1) 2 -B r = lME r . At t = 1ms, 
the soliton energy has decreased already to Q.96E r , the 
rest of energy has been taken by the expelled atoms. A 
steady state is reached for a soliton energy of E = 0.92EV. 
While moving, some energy is devoted to cross the PN 
barrier (the soliton configuration changes its shape from 
a configuration centered in one well of the optical lattice 
to a configuration centered in one maximum and vicev- 
ersa) . This change of the shape of the soliton is reflected 
in the out of phase oscillations of the kinetic energy with 
respect to the potential plus non linear energy in such a 
way that the mean value of the energy remains constant. 



III. "EFFECTIVE" POTENTIAL BARRIER 

We discuss first the interaction of a bright lattice soli- 
ton with an "effective" potential barrier. Scattering 
depends on the width of the defect (/) and the rele- 
vant energy scale, settled by the ratio \V m \/Ek, where 
Ek =< P > 2 /2m is the fraction of the total kinetic en- 
ergy Ej. devoted to move the soliton, and <> denotes 
time average (before reaching the defect). The momen- 
tum P is defined as: 



P(t) 



-itup*(x, t)\7ip(x, t)dx. 



(6) 



The rest of the kinetic energy is needed to keep the struc- 
ture and cannot be used to overcome the "effective" po- 
tential barrier. We have checked that apart from the 
necessary change in shape to overcome the PN barrier, 
the soliton keeps its overall shape when it reaches the 
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FIG. 1: (a) Transmission coefficient, T, as a function of the 
amplitude of the "effective" potential barrier, \V m \; (b) Trans- 
mission coefficient, T and (c) reflection coefficient R as a func- 
tion of time for different amplitude defects: \V m \ = (solid 
line), \V m \ = 0.01£; r (dashed line), \V m \ = 0.0115B,- (dot- 
dashed line), \V m \ = 0.0117-E r (circles) and \V m \ = 0.018£ r 
(dotted line). In all the plots the soliton kinetic energy is 
E k = 0.01£V and the width of the defect I = 2d 



defect. This corroborates that there is no transfer be- 
tween non linear energy and kinetic energy apart of the 
one corresponding to the already discussed PN barrier. 

We distinguish two regimes of parameters: (i) when 
the amplitude of the "effective" potential barrier is on 
the order of the kinetic energy of the soliton (\V m \ ~ Ek) 
and (ii) when amplitude of the potential barrier is much 
larger than the kinetic energy of the soliton. In the for- 
mer case, the potential barrier acts as a quantum switch, 
i.e., either the entire soliton is transmitted or it is com- 
pletely reflected depending on the amplitude of the bar- 
rier (Fig. [3Ja)). The transmission (T) and reflection (R) 
coefficients are calculated by integrating over space (and 
time) the density of the wavefunction in the region after 
and before the defect, respectively. Note that since only 
approximately 35% of the initial atoms form a soliton and 
since there are also loses of atoms during the kick, the 
merit figure for perfect transmission is well below 1 and 
corresponds approximately to 0.27 (N=135). For a fixed 
width of the defect a drastic change of behavior occurs for 
a given height of the barrier \V£\. The wider the defect 
is the lower the critical value of the height of the barrier 
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FIG. 2: (a) Reflection and transmission behavior of a soliton 
interacting with an effective potential barrier as a function of 
the potential width I (in d units) and potential amplitude \ V m 
(in E r units): the region of total transmission in gray, and the 
region of total reflection in white. Inside the reflection region, 
bands in which trapping and reflection occurs appear. Region 
of parameters for which tunneling and overbarrier reflection 
occur are also shown. The dashed line shows the value of 
the kinetic energy of the soliton. (b) Density profiles of the 
trapped structure that appears for \ V m \ = 0.55-E r (lower plot), 
\V m \ = 0.8E r (middle plot) and \V m \ = l.l-E r (upper plot). 
In the three cases I = 8d. 



\V^\. The critical values, indicating the transition be- 
tween complete transmission and complete reflection for 
different potential widths are shown in FigEJa) by solid 
black squares. Below these values, depicted by a gray 
region in FigOJa), complete transmission of the soliton 
occurs. As one approaches the critical value from below, 
the soliton experiences a time delay with respect to free 
propagation (i.e., in the absence of the defect). This de- 
lay increases as one gets closer to the critical point and 
eventually the time needed by the soliton to cross the 
barrier diverges (see Fig. [3]). Above the transition line, 
reflection of the entire soliton occurs after a storage time 
inside the region of the barrier that also increases as one 
approaches the critical value. To illustrate this behavior, 
Fig. [ljb,c) shows transmission and reflection coefficients 
as a function of time for a barrier of fixed width I = 2c? 
and different values of the amplitude. In the situation 
shown in Fig. [TJ the critical value is nearly equal to the 
kinetic energy of the soliton but if the width of the bar- 
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FIG. 3: Transmission delay time as a function of \V m \ for a 
soliton with kinetic energy Ek = 0.01E r crossing a defect of 
width I — d 



rier is reduced, this critical value can exceed the kinetic 
energy of the soliton. In this case, the soliton tunnels 
through the barrier, i.e., transmission is obtained for val- 
ues of the amplitude of the barrier higher than the kinetic 
energy of the soliton (see Fig. Eta)). On the other hand, 
for wider defects, a region of overbarrier reflection ap- 
pears (Fig. [UJa)). There, the lattice soliton is completely 
reflected although it has a kinetic energy larger than the 
height of the potential barrier. This region extends for 
a wide range of widths of the defect. We have checked 
that overbarrier reflection occurs even in the limit when 
the width of the defect is much larger than the size of the 
soliton. 

Fig. [3] shows the delay time in transmission with re- 
spect to the absence of defect, td, for I = d as a function of 
the height of the barrier, \V m \, including the cases where 
tunneling occurs (\V m \ > 0.01E r ). 

Up to now we have described the transition from com- 
plete transmission to complete reflection by fixing the 
width of the defect and varying the amplitude. It is 
worth noticing that a similar switching behavior can be 
obtained by fixing the amplitude of the effective poten- 
tial and changing its width. This would correspond to 
horizontal lines in Fig. [2fa) crossing the transition line 
(solid black squares) in the region where \V m \ is on the 
order of the kinetic energy of the soliton. This observed 
abrupt transition from complete reflection to complete 
transmission opens the possibility to use the system as a 
quantum switch. A similar switching behaviour has been 
predicted for optical Bragg solitons described with the 
coupled mode equations [18j |. 

Let us now focus on the regime where the amplitude 
of the barrier is much larger than the kinetic energy of 
the soliton \V m \ ^S> Ek, where the expected behavior is 
complete reflection of the soliton. Complete reflection 
occurs but there are specific values of the ratio i/|K„| for 
which the soliton splits into two parts: a fraction of the 
initial soliton becomes trapped inside the region of the 
barrier while the other part is reflected back keeping a 
solitonic structure. Fig. [U^a), for \V m \ >• Ek, shows the 
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regions where the soliton splits into two parts (trapping 
and reflection) embedded in the complete reflection re- 
gions. The fraction of atoms trapped inside the defect 
region has its origin on the atoms lost by radiation due 
to the repulsive force experienced by the soliton when it 
reaches the potential barrier. These radiated atoms enter 
the region of the barrier and for some specific ratios of the 
width and the height of the defect the fraction of trapped 
atoms increases. These trapping regions appear as bands 
as shown in Fig. Ufa). In each band, the trapped frac- 
tion exhibits different spatial distributions: for the first 
(lowest) one, the structure is a single hump; in the second 
one a double hump structure appears, and so on (see Fig. 
Efb)) . A noticeable feature of this trapped fraction is that 
the density maxima of the structure are located at the po- 
sitions of the maxima of the optical potential. Increasing 
the amplitude of the barrier, the structure becomes more 
independent of the lattice periodicity. The extension of 
the trapped structure is the same independently of the 
features of the barrier but the number of trapped atoms 
differs for different widths of the barrier. The narrower 
the defect is the larger the number of trapped atoms. 
This number changes also with \V m \ inside each band, 
being maximum at the center of the band. For all cases 
the number of atoms forming the reflected soliton is al- 
ways larger than the trapped fraction. We have checked 
that these "resonance" bands like do not correspond to 
bound states of the linear case. We have also observed 
that this behaviour occurs for all the accessible initial 
transfers of momentum that allow motion of the soliton. 



IV. "EFFECTIVE" POTENTIAL WELL 

Let us now turn to the interaction of a lattice soli- 
ton with an "effective" potential well with a depth of 
the order of its kinetic energy. For a fixed depth of the 
well, the soliton exhibits different behaviors depending 
on its kinetic energy. For low kinetic energies, the soli- 
ton gets bound with the defect and exhibits oscillations 
while for kinetic energies overcoming a certain thresh- 
old, the soliton crosses the defect region. In the latter 
case, the only detectable effect of the potential well is the 
speed up of the soliton with respect to free propagation. 
It is important to note that as the width of the defect 
increases, the range of velocities for which transmission 
occurs decreases. To illustrate the described behavior, 
Fig. HJa) shows the time evolution of the trapped frac- 
tion density, D, for different initial transfer of momen- 
tum: p = 0.05kh (solid line), p — Q.lkh (dashed line), 
p = O.I7fcft (dotted line) and p = 0.2kh (circles) keeping 
the depth (\V m \ = 0.018EV) and the width (I = 8d) of 
the well fixed. The positions in time of the minima of 
the trapped fraction correspond to the turning points of 
the oscillating movement of the soliton around the "ef- 
fective" well. The maxima indicates the times for which 
the soliton is completely inside the well. As expected, the 
amplitude of the oscillations increases with an increasing 
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FIG. 4: (a) Temporal evolution of the trapped fraction of the 
soliton interacting with an "effective" well of depth \V m \ = 
0.018i5 r and width I = 8d after an instantaneous transfer of 
momentum of p — 0.05kh (solid line), p = Q.lkh (dashed line), 
p = 0.17kfi (dotted line) and p — 0.2kh (circles), (b) Contour 
plot of the evolution in space and time of the lattice soliton 
with conditions corresponding to the dashed line case in (a). 



momentum transfer. If the amplitude of the oscillation 
is larger than the width of the defect, the turning points 
are located outside the potential well. This is reflected 
by a lower value of D. Fig. 01b) shows a contour plot of 
the evolution in space and time of a lattice soliton with 
E k = 0.01E r (dashed line case in Fig. Ufa)). The width 
of the " effective" well is shown at the right hand part of 
the plot to illustrate that indeed the turning points are 
outside the defect. 

By fixing Ek = 0.01£V, we explore now the depen- 
dence of the oscillations on I and \V m \. Fig 03(a) displays 
the temporal evolution of the trapped fraction, D, for 
\V m \ = 0.018i? r and different values of the width of the 
defect: I — Ad (solid line), I = 8d (dashed line) and 
I = 12d (dotted line). The frequency of the oscillations 
gives an indication of the width of the defect, decreasing 
as the width increases, while the amplitude remains ap- 
proximately the same for all widths. In Fig. [5Jb) we fix 
I = 4d and display the trapped fraction, D, as a function 
of time for \V m \ = 0.018£V (solid line), \V m \ = 0.03E r 
(dashed line) and \V m \ — 0.08_E r (dotted line). By in- 
spection of Fig. E^b), one can confirm that the frequency 
of the oscillations increases with the depth of the po- 
tential while the amplitude of the oscillations decreases. 
This is due to the fact that the soliton experiences a much 
larger attractive force as the depth of the defect increases 
limiting the displacements around the central position of 
the well. 
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FIG. 5: Temporal evolution of the trapped fraction density 
of a lattice soliton with = 0.01-E r interacting with an 
"effective" well with (a)|V m | = 0.018£ r and / = Ad (solid 
line), I = 8d (dashed line) and I — 12d (dotted line); (b) I = 4<f 
and \V m \ = 0.018£ r (solid line), \V m \ = 0.03£ r (dashed line) 
and \Vm\ = 0.08E r (dotted line) 



The trapping of the entire lattice soliton around the po- 
sition of the defect opens possibilities to use the system 
as a quantum memory because it provides the capacity of 
storage. Nevertheless, in order to perform a memory, one 
should also be able to release the trapped structure after 
a desirable time and with the minimum loses. We have 
checked that a soliton trapped in an "effective" poten- 
tial well can be released with a certain velocity keeping 
the totality of its initial atoms if the defect amplitude 
is instantaneously set to zero. In fact the velocity of 
the lattice soliton after releasing it will depend on the 
amplitude of the oscillations it was performing while it 
was trapped. Specifically, the velocity of the structure, 
after releasing it, grows with the amplitude of the os- 
cillations. Moreover, choosing appropriately the time at 
which the release takes place, one can vary the direction 
of the movement. 



V. CONTROL OF THE COLLISIONS 

Now we investigate if the inclusion of a defect in the lat- 
tice helps to control the interactions between two lattice 
solitons. It has been shown that collision of two identi- 
cal lattice solitons (moving with the same velocity and 
with the same average phase) merge into a soliton with 
the same number of atoms as the initial ones [9]. The 
exceeding atoms are lost by radiation. If an "effective" 




20 40 60 80 100120140 

t(ms) 



20 40 60 80 100120140 

t(ms) 



20 40 60 80 100120140 

t(ms) 



FIG. 6: Contour plot of the evolution in space and time of 
a collision between two identical solitons initially placed at 
symmetrical positions with respect x — and moving in op- 
posite directions when an "effective" potential barrier with 
I = 2d and (a) \V m \ = 0.012S r , (b) \V m \ = 0.2E r and (c) 
\Vm \ = 0.5E r is placed at x = 0. 



potential barrier much narrower than the dimensions of 
the solitons is placed at the crossing point, we find the 
following behaviours: (i) for \V m \ < the merging be- 
havior is maintained (Fig. [5] (a)); (ii) for \V m \ » 
each soliton reflects back (Fig. [6] (b)). Moreover, for 
some values of \V m \, in addition to the reflection, a frac- 
tion of atoms trapped in the defect can appear (Fig. [5] 
(c)). The trapped fraction shows up the same features as 
in the single soliton case (section III). Modifying the fea- 
tures of the defect, different outcomes can be engineered. 
For instance, when the width of the barrier is of the order 
of the dimensions of the initial solitons, effects like the 
trapping of both solitons at the edges of the barrier can 
occur. 



VI. CONCLUSIONS 

Summarizing, we have found that bright matter wave 
lattice solitons behave as "quantum" particles when col- 
liding with an "effective" barrier/well, corresponding to 
a defect in the optical lattice. Among the rich dynamics 
exhibited by the system, we would like to remark two 
effects. The first one corresponds to the interaction of 
a soliton with an " effective" potential barrier which per- 
mits the implementation of a quantum switch. In this 
case, a sharp transition from complete reflection to com- 
plete transmission is present at a specific value of the 
height of the barrier. Although this resembles the classi- 
cal particle behaviour, the quantum nature of the solitons 
is explicitly manifested in the appearance of overbarrier 
reflection and tunneling. The second effect we would like 
to stress appears when the defect acts as an "effective" 
potential well. We have shown that trapping of the entire 
soliton around the position of the defect and its release 
on demand with a given velocity and direction of mo- 
tion is possible. This fact indicates the suitability of the 
system as a quantum memory. Finally, it has been also 
reported that the presence of a defect in the lattice can 
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help to control the interactions of two lattice solitons. 
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